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Abstract High-quality cross sections for the reaction — > ir d have been measured using the CLAS 
at Jefferson Lab over a wide energy range near and above the 77-meson photoproduction threshold. At 
backward cm. angles for the outgoing pions, we observe a resonance-like structure near E 7 =700 MeV. Our 
model analysis shows that it can be explained by r\ excitation in the intermediate state. The effect is the 
result of the contribution of the iV(1535)Sn resonance to the amplitudes of the subprocesses occurring 
between the two nucleons and of a two-step process in which the excitation of an intermediate r\ meson 
dominates. 



PACS. 13.60.Le Meson Production - 14.20.Gk Baryon resonances with S = - 25.20.Lj Photoproduction 
reactions 
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1 Introduction 

Interactions of the 77 meson with few-nucleon systems com- 
plement our knowledge of the 77-nucleon interaction. Inter- 
est in these systems stems from the hypothetical existence 
of ry-nuclear quasibound states. Such states have been pre- 
dicted by Haider and Liu [1] and Li, Cheung, and Kuo [2] . 
Although there has been no direct experimental verifica- 
tion of this hypothesis to date, there is mounting evidence 
that such states might exist in the lightest few-nucleon 
systems [3-5]. 

For the case of the three-nucleon system, it was found 
at Saclay that the dp — * r\ 3 He production amplitude falls 
rapidly just above the r\ threshold [3]. A less pronounced 
slope was found in the dd — > r\ 4 He amplitude [4] . For the 
two-nucleon system, very strong final-state interactions 



(FSI) were found in the pp —> pprj cross section in the 
threshold region [5] . The energy dependence of the NN — > 
NNrj and NN — ► drj reactions can be understood in terms 
of NN FSI ( [6] and references therein) . However, as noted 
in Ref. [7], the existence of a narrow virtual state in the 
ry-deuteron system can be inferred. 

The production of a virtual r\ meson may also play a 
role in other nuclear reactions, even in those for which 
there is no r\ in either the initial or the final state, but 
only in an intermediate state. Examples are the reactions 
pd — ► 7T + 3 H and pd — > 7r° 3 Hc, for which there are strong 
indications that an intermediate r\ cusp is present [8] . Ev- 
idence for an intermediate r\ meson also was found in elas- 
tic ird backward scattering. This contribution manifests 
itself as a cusp in the energy dependence of the backward 
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differential cross sections near the rj threshold. The effect 
was predicted in Ref. [9] and was confirmed by several in- 
dependent measurements of backward ird scattering [10]. 
In this work, we present the first systematic evidence for 
a similar phenomenon (the first indication was found in 
Ref. [11]) in the coherent photoproduction of a neutral 
pion on the deuteron, 7c?— >ir°d. 

All of the above phenomena take place because r\- 
meson production in hadron-hadron collisions near thresh- 
old is enhanced. This is because the cross section for ex- 
citation of the nearby baryonic resonance 7V(1535)Sii is 
large and this resonance is strongly coupled to the r\N 
channel. Since the amplitudes for photoproduction of the 
iV(1535)5n are also large [12-14], one can expect a similar 
enhancement in various photonuclear reactions. 

Recently, coherent photoproduction of the 7r° meson 
from the deuteron was studied theoretically [15]. In par- 
ticular, it was demonstrated that at large cm. scatter- 
ing angles and photon energies E 1 between 600 and 800 
MeV, the two-step process with the excitation of an in- 
termediate 77-meson (shown in Fig. H^b)) dominates over 
single-step photoproduction (shown in Fig.QTa)) and pion 
rescattering. This two-step process can be analyzed as two 
sequential subprocesses: 7-ZV1 — > 7V(1535) — * r\N\ and 
r\N 2 -> V(1535) -> 7r iV2, where N x and N 2 are the 
two nucleons in the deuteron. It was shown in [15] that 
this mechanism explains qualitatively the structure in the 
~fd^ir°d differential cross section, which we present here, 
at large angles and for E 7 ~ 600-800 MeV. The main 
conclusions of Ref. [15] were reproduced in a more recent 
paper [16], where it was shown that in addition to this 
two-step process, the full dynamics in the intermediate 
NNr] system could be important as well. Other theoret- 




a) M a h) M b 



Fig. 1. Feynman diagrams for the yd—>ir d reaction consid- 
ered in [15]: (a) single-scattering amplitude M a ; (b) double- 
scattering amplitude Mi,. It was shown in [15] that (b) domi- 
nates over (a) at backward angles for E~, ~ 700 MeV. 

ical studies of the yd — ► ir°d reaction can be found in a 
number of papers [17]. However, none of these considers 
the effect of the opening of the r\ threshold at 700 MeV. 

Our photoproduction data, presented here, give for the 
first time clear evidence for a prominent effect around 700 
MeV at large cm. angles, which can be explained by the 
excitation of an intermediate 77 meson. 



2 Experiment 

The experiment took place in Hall B at Jefferson Lab using 
the CEBAF Large Acceptance Spectrometer (CLAS) [18]. 



A collimated, tagged, real-photon beam, with energies be- 
tween 0.5 and 2.3 GeV was incident on a 10-cm-long LD2 
target installed in the center of the CLAS. The energies 
of the photons were tagged using the Hall-B tagger [19]. 
The outgoing deuterons were tracked in the six toroidal 
magnetic spectrometers of the CLAS. They were bent out- 
wards by the magnetic field, and their trajectories were 
measured by three layers of drift chambers surrounding 
the LD2 target. The time of flight of the deuterons was 
measured by 6x48 scintillators (TOF) that surround the 
CLAS detector outside of the magnetic field. A set of 
six scintillator counters, comprising the start counter and 
placed just around the target, measured the event time at 
the vertex. The CLAS covers the polar angular range be- 
tween 8° and 142° in the laboratory system and the entire 
range in azimuthal angle. 



3 Data Analysis 

Since the CLAS has very good acceptance for charged 
particles, and a very limited one for neutrals, the anal- 
ysis of the 7c? — > ir°d reaction was based on detecting 
the final-state deuterons and selecting the good events 
by the missing-mass technique. Deuterons were identified 
from their time of flight and momentum, which allowed a 
mass reconstruction. Thus, the initial event selection was 
done based on the reconstructed mass. Further selection 
of events was achieved by comparing the event vertex time 
with the photon vertex time as measured by the tagger. 
Low-momentum protons were discarded based on a cut 
on the particle momentum vs. energy loss in the TOF. In 
addition, fiducial cuts were applied to the remaining data 
sample in order to remove edge areas of the detector where 
the acceptance was not well reproduced by a simulation. 

Once the data were reduced, based on all of the above 
cuts, they were binned in photon energy and pion cm. 
scattering angle (E 7 , cos 9* ) . Here, we present differen- 
tial cross sections for the reaction yd —> ir°d based on 
a photon-energy bin width of 25 MeV and a cos 8* bin 
width of 0.1. For every such (E^,cos6*) bin, we deter- 
mined the reaction yields after sideband background sub- 
traction was performed on the missing-mass (mrnj) distri- 
butions as shown in Fig. [5J The quantity mm 2 , is defined 
as mm? d — (p 7 + Pt ~ Pd) 2 , where p 7 , pt, and pd are the 
four-momentum vectors of the beam, the target, and the 
recoil deuteron, respectively. The systematic uncertainty 
associated with the background subtraction depends on 
the kinematic bin, and varies between 1.5% and 6.4%. 
The statistical uncertainty of the extracted yields is also 
bin-dependent, and varies between 3% and 30%. 



4 Results 

Differential cross sections were obtained by normalizing 
the true-event yields to the photon flux, the number of 
target scattering centers, and the CLAS acceptance [20]. 
The statistical uncertainty of the photon flux is negligible. 
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Fig. 2. An example of the background subtraction, for the bin 
0.675 < E 1 < 0.7 GeV and -0.8 < cos9* < -0.7. One sees 
that the pion peak is well separated from the background. The 
contribution of the latter to the peak is < 10%. The shaded 
areas show the sidebands around the pion peak we used in order 
to determine the shape and the magnitude of the background 
underneath it. 



The systematic uncertainty of the evaluated photon flux 
is 3.3%. The statistical uncertainty associated with the 
calculated CLAS acceptance is bin-dependent and varies 
between 0.8% and 2.5%. The systematic uncertainty of 
the acceptance is less than 10%. A common factor of 
1.0141 ± 0.0006 stat ± 0.0005 syst , due to the inefficiency 
of the procedure for choosing the right photon for a given 
event, was applied to all of the differential cross sections. 
There is also an overall systematic uncertainty of 0.5% 
related to the determination of the target length and den- 
sity [21]. Thus, the total uncertainty of the differential 
cross sections is bin-dependent and varies smoothly from 
11% to 33%. 

Our differential cross sections for the reaction 'yd-^TT°d 
for five cm. angular bins are shown in Fig. [3l Note that 
some of the angular bins overlap partially with each other. 
The latter is due to the fact that for a consistency check, 
wc determined differential cross sections for two entirely 
separate angular binnings. The figure illustrates the con- 
sistency of the structure and the model interpretation at 
different kinematical binnings. 

Overall, our data are in very good agreement with 
previously measured differential cross sections [11,22-24]. 
However, in the range of photon energies where we observe 
the backward peak discussed here, the data of Ref. [11] ex- 
hibit a structure but of a smaller magnitude, whereas the 
data of Ref. [23] do not exhibit any structure. In order to 
understand this discrepancy we performed many consis- 
tency checks, and our analysis procedure was studied in 
great detail for possible sources of errors [20]. The struc- 
ture near 0.7 GeV and its magnitude persists. Both of 
the measurements [11,23] were done at the same acceler- 
ator facility and used an untagged bremsstrahlung beam, 
which might introduce large systematic uncertainties in 
the determination of the photon flux. 




0.5 0.6 0.7 0.8 0.9 
E y (GeV) 

Fig. 3. Excitation functions for the reaction "fd — > ir°d for sev- 
eral bins in cos 9* . The filled circles are our CLAS data. The 
error bars show the total bin-dependent uncertainty for each 
data point. The dashed and solid curves are the results of fit 1 
and fit 2 (see the text), respectively, obtained with the helicity 
amplitudes from Ref. [13] (set (2)). One sees that the excita- 
tion function at cos 9^ = 0.0 does not exhibit any prominent 
structure near 0.7 GeV. Note that some of the angular bins 
overlap partially with each other. 



5 Discussion 

To achieve a quantitative understanding of our data, we 
employ the semiphenomenological description of Ref. [15]. 
We express the reaction amplitude M as a sum M = 
MiS +M 2 S, where Mi S is the two-step amplitude M b 
given by the diagram in Fig.QJb) and calculated in Ref. [15] 
with intermediate r\ production, and S is a spin factor. 
M2S is the effective "background" amplitude which takes 
into account all possible background diagrams (including 
the single-scattering amplitude M a shown in Fig. [31 a)). 
We parametrize it as M 2 = Aexp(i(f2 — bE 7 ), where if2 is 



5 



the phase. The square of the total amplitude with unpolar- 
ized particles is then written as |M| 2 = \M 1 +M 2 \ 2 \S\ 2 and 
is applied to describe the experimental excitation func- 
tions for given values of cos 9* . We give the results of fits 
for two different parametrizations of if 2' 

fit 1 : ^2(^7) = ol + /3(E 7 -0.7 GeV); (1) 
fit 2: ip 2 {E 1 ) = a + <pi(Ey). (2) 

For fit 1 , we use a linear independent background phase 
(f2(E 7 ) with two parameters, a and (3. We define a = 
Lp2(Ej) at energy E 1 — 0.7 GeV, where "by eye" the am- 
plitude M\ peaks. Thus, we use four parameters (A, b, a, 
and 0) in fit 1. 

For fit 2, we use the parameter a as the relative phase 
of the amplitudes M\ and M%, i.e., a = tp2{Ey) — ipi(E 7 ) , 
where (pi (E 1 ) is the phase of the amplitude Mi and is a 
given function. Thus, we use three parameters (A, b, and 
a) in fit 2. In both variants all the parameters were varied 
independently to fit the data for various values of cos 9*. 

We note that the two-step amplitude Mi with inter- 
mediate r\ production is proportional to A^ 2 — A™^ 2 , where 

^1/2 (^1/2) 1S ^ ne helicity amplitude of the decay iV(1535) - 
£7(777). Values for these amplitudes vary widely in the lit- 
erature. Here, we give some sets (in units of 10 -3 GeV -1 / 2 ): 

(1) A\ l2 = 107, A» /a = -96 ( [14]); 

(2) A\ j2 = 78, A™ /2 = -50 ( [13], used in [15]); 

(3) A\ j2 = 60, A n lj2 = -20 ( [25]). 

The results for the free parameters from the data fit- 
ting are given in Table [1] for large values of 9*. Here, the 
coupling set (2) was used. The dashed and the solid curves 
in Fig. 3 correspond to the results of fits 1 and 2, respec- 
tively, and both fits give a satisfactory description of the 
data. The corresponding values of x 2 /iV (N is the num- 
ber of degrees of freedom) are also given in Table [TJ Using 
this procedure, we obtain a good description of the data at 
large scattering angles (cos 9* < —0.5), where the r\ effect 
is strongly pronounced. We do not consider the data at 
smaller angles because a more complicated parametriza- 
tion of the background is needed and in any case the effect 
is less pronounced. 

The backward-angle structure is qualitatively repro- 
duced using various sets of couplings A^ 2 and A^ 2 . We 
obtained also a good description (not shown) using cou- 
pling set (1) in fit 2 and set (3) in fits 1 and 2. Thus, 
the experimentally observed enhancement in the excita- 
tion functions near E 1 ~ 0.7 GeV needs no exotic expla- 
nation such as was done, for instance, in Ref. [11]. 

In order to evaluate the statistical significance of the 
structure exhibited by our data, we performed a x 2 test of 
goodness of fit. We tested the hypothesis that the excita- 
tion functions at large angles, which are shown in Fig. [3l 
can be described by a smooth function. We used the am- 
plitude M2, which falls off smoothly with E 1 and describes 
well the experimental data outside the energy range where 
the enhancement is observed, to fit the excitation func- 
tions at large angles. The results for the x 2 an d the signif- 
icance level of this hypothesis (the P-value) are given in 



Table 2. \ 2 an d number of degrees of freedom N of the fits to 
the large angle excitation functions with the smooth amplitude 
Mi. The significance of the hypothesis (the P-value for the \ 2 
test of goodness of fit) that M2 describes the data is also shown. 



cos 8* 


x 2 


N 


P-value 


-0.7 


44.43 


22 


0.003 


-0.75 


43.68 


22 


0.004 


-0.8 


35.71 


21 


0.02 


-0.85 


77.07 


22 


0.0001 



TableO We see that the probability that this hypothesis is 
true is smaller than 2%. Thus, the statistical significance 
of the structure is greater than 98%. 

Finally, the excitation of the 7V(1535)Sii resonance de- 
pends on the initial polarization of the beam and target. 
The spin factor S of the two-step amplitude with an inter- 
mediate rj meson and an S'-wave deuteron wave function 
is S = (£2 ' [ e x where e and ei (£2) are the polariza- 
tion 3-vectors of the initial photon and the initial (final) 
deuteron. For an unpolarized final-state deuteron and for 
A 7 = ±1 and A^ = 0, ±1 as the helicities of the initial 
photon and deuteron, respectively, we obtain \S\ 2 = 2/3 
for unpolarized particles; |5| 2 = 1 for = 0, — A 7 ; and 
\S\ 2 = for Arf = A 7 . Thus, for A^ = 0, — A 7 , we expect 
the 77-effect to be enhanced compared with the unpolar- 
ized case. For A^ = A 7 (parallel polarization of the initial 
photon and deuteron), the excitation of the ^(1535)511, 
and hence the 77-effect, should be suppressed. 



6 Summary 

To summarize, we have measured unpolarized differential 
cross sections for the -fd^ir°d reaction at backward cm. 
angles and for photon energies above 500 MeV. The data 
show a pronounced structure (with statistical significance 
greater than 98%) in the excitation functions in the re- 
gion of 700 MeV. For the first time, this phenomenon is 
systematically studied with good accuracy. The structure 
can be explained by the opening of the ?7-photoproduction 
threshold on a single nucleon, and in our model is mainly 
related to the excitation of the intermediate resonance 
V(1535)Sii in a two-step process. The details of the un- 
derlying dynamics can be further explored via polarization 
measurements. 
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Table 1. Parameters and \ 2 P er degree of freedom N of the model-based fits, explained in the text, to the experimental 
excitation functions at large cm. angles. Parameters A, b, a, and /3 (A, b, and a) were used in fit 1 (fit 2). The values for the 
helicity amplitudes A^'Z used in the calculations were taken from Ref. [12] (set (2)). N = Ndat — N par , where N par = 4(3) for 
fit 1(2). N dat = 24 for cos 6»* = —0.8, and 25 for the others. 



cos 0% 


fit 


A 


b (GeV" 1 ) 


a (deg) 


(3 (deg/MeV) 




-0.7 


1 


1.68 


2.53 


177 


0.79 


1.12 


-0.7 


2 


1.55 


2.40 


59 




1.73 


-0.75 


1 


1.57 


2.51 
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0.80 


1.27 


-0.75 


2 


1.52 


2.42 


59 




1.07 


-0.8 


1 


1.69 


2.68 


158 


0.81 


1.55 


-0.8 


2 


1.63 


2.55 


66 




0.83 


-0.85 


1 


1.51 


2.59 


136 


0.89 


0.99 


-0.85 


2 


1.55 


2.61 


42 




0.79 
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